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EFFECTS OF MISTUNINO ON BENDING-TORSION FLUTTER AND RESPONSE OF A CASCADE TN INCOMPRESSIBLE FLOW 

by 

Krishna Rao V» Kaza^ 

Th^ University oi Toledo and 

National Aeronautics and Space Administration 
Lewis Research Center 
Clevelondi Ohio AA135 

end 

Robert E. Klolb** 

National Aeronautics and Space Administration 
Lewis Research Center 
Clevel^’nd^ Ohio AA13B 


Abstract 

This paper presents an Invoetigatlon of the 
effects of blade mistuning on the oeroelostic sta- 
bility end response of a cascade in incompressible 
flow* The aerodynamic^ inertial, and structural 
coupling between the bending and torsional motions 
of each blade and the aerodynamic coupling between 
the blades are included in the formulotion* A dig- 
ital computer program was developed to conduct 
parametric studies* Results Indicate that the mis- 
tuning has a beneficial effect on the coupled 
bending- torsion and uncoupled torsion flutter. The 
effect of mistuning on forced response, howeve’', 
may be either beneficial or adverse, dopendinL on 
the engine order of the forcing function. Addi- 
tionally, the results illustrate that it may be 
feasible to utilize mistuning as a passive control 
to increase flutter speed while maintaining forced 
response at an acceptable level. 


Nomenclature 

[A] aerodynamic matrix due to motion 

[Ar] aerodynamic matrix due to motion In 

rch mode; r » 0, 1, 2 ... N-1 

{AD} aerodynamic matrix due to wake induced 

flow 

{AD,^} aerodynamic matrix due to wake induced 

flow in the rth mode, r « 0, 1, 2 ... 
N-1 

a elastic axis location, nondimensional 

b semichord 

c chord 

[D] ,[Dq] matrices defined in equation (11); 

s « 0, 1, 2 ... N-1 

[E] matrix defined In equation (4) 

E(s,r) defined in equation (A) 
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tG],[Ge] 

%hg-CKaa 

ha 

har 

[I] 


1 

ha O.B 


kp 

tM 




^hhr»^har 

^ahr»^aar 

^whr>^war 


base for noturak logarithm 

matrices defined in equation (11); 
s ^ 0, 1, 2 ... N-1 

quantities defined in equation (11) 

b*^nding deflection of sth blade 

bending deflection of blade In rth mode 
of tuned cascade 

unit matrix 

mass moment of inertia of sth blade 
about elastic axis per unit spun; 


V=r 


’^Adjunct Professor, Mechanical Engineering Depart- 
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bending and torsional stiffness respec- 
tively, of sth blade 

reduced frequency, uib/V 

reduced flutter frequency, wpb/Vp 

lift due to motion of sth blade per 
unit span, positive up 

lift due to wakes of sth blade per unit 
span, positive up 

nondimensional lift coefficients due to 
bending and torsional motions, respec- 
tively, in rth mode 

nondimensional moment coefficient due 
to bending and torsional motion, re- 
spectively, in rth mode 

nondimensional lift and moment coeffi- 
cients, respectively, due to wake in 
rth mode 

moment about the elastic axis due to 
motion of sth blade per unit span, 
positive nose up 

moment of sth blade per unit span about 
the elastic axis due to wake, posi- 
tive nose up 
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rOg maos per unit upon pf oth blade 

N number of bladeo in caoeode 

[PJ matrix defined in equation (ll) 

*^hcr*'^hor nmXtibXado coordlnateo for bending 
motion in rth mode 

mulciblado coordlnateo for toroional 
motion in rtlt mode 

r integer specifying the mode of tuned 

rotor; r « 0# li 2 ♦.* N“l| also the 
engine order of the excitation 

tg radius of gyration of sch blade, non- 

® dlmonsionallaed with respect to b 

a integer specifying blade, a « 0, I, 

2 ... N-l; also blade spacing (fig* 1) 

Sa static mass moment of ach blade per 

® unit span about elastic axis, posi- 

tive when center of gravity la aft of 
elastic axis 

t time 

V frcGstream velocity relative to the 

blade 

Vj? flutter speed 

Wj- velocity induced by wakes 

(X) column matrix, defined in equation (4) 

X,Z rectangular coordinate axes 

Xq. dimensionless static unbalance of sth 

® blade (“Sd^/rngb) 

{Y} column motrlx, defined in equation (4) 

ttg amplitude of torsional motion of sth 

blade, positive clockwise 

as, Id torsional amplitude of each blade of 

tuned rotor 

amplitude of torsional deflection of a 
blade in rth mode of a tuned cascade 

Sr interblade phase angle, 2irr/N 

p 

Y nondimensional eigenvalue, (wg/w) 

Yhg nondimensional uncoupled bending fre- 

quency of sth blade 

Yojg nondimensional uncoupled torsional fre- 

quency of sth blade 

6hg>‘Sag logarithmic decrements of sth blade in 

bending and torsion, respectively 

Ch >^0 damping ratios of sth blade in bending 

® ® and torsion, respectively 

n location of elastic axis measured from 

leading edge, (a + l)/2 



a^Jmuthal position of sth bXado» do^ 
fined in equation (6) 

H 

masp do of ath blade, 

V 

real part of eigenvalue, defined in 
equation (12) 

V 

imaginary part of eigenvalue, defined 
in equation (12) 


nondimensional flutter frequency 


stagger angle, figure I 

P 

fluid density 


frequency 

Wg 

referenco frequency 

“•'8 

yKh^/mg 

“ag 


1 1.( } 

matrices 

(') 

differentiation with time 

1 ]"^ 

Inverse of a matrix 

E 

indicate summotion over r 0, 1, 2 . 
N-1 


I. Introduction 

In the development of modern aircraft turbofan 
engines, the aoroelastlc stability and response of 
bladed-disk assemblies have been among the most 
difficult problems encountered. The study of sta- 
bility and response in these assemblies is compli- 
cated by the presence of small differences between 
the Individual blades, known os mistuning# The 
published results in this area which will be dis- 
cussed later have shown that mistuning can have a 
beneficial effect on turbine engine blade flutter 
and on adverse effect on forced response. Experi- 
enced* 2) has further shown that there have been 
costly failures in the development and production 
phases in which mistuning appeared to have played 
an Important role. Thus, an improved basic under- 
standing of these effects is important in the 
design phase. 

To improve the basic understanding of the ef- 
fects of mistuning on aeroelastic stability and 
response and then to explore the possibility of 
utilizing mistuning as a passive control to allev- 
iate flutter and to minimize forced response, an 
effort has been in progress in the Structural 
Dynamics Section, NASA Lewis Research Center. As a 
part of this general effort, the effects of mis- 
tuning on coupled bending-torsion flutter and on 
aeroelastic response due to wakes have been studied. 
This paper presents the results of the study for 
Incompressible flow. 

Either because of the complexities or because 
of the general belief that the turbomachinery blade 
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flufcccr Involvoo only a single dogroo of froo*- 
dotn, previous reoenreherot^"^) have studied the 
effect of blade miotuntng on aoroolaoclc atobility 
mid rooponao by conoiderins either pure bending 
ttKxClon or pure torolonal mociem of the blmleo. The 
itudloQ of mi itl-dogroo of freedom blade flutter In 
vbe publiohed lltorature(5»8»9»l0ill) have been 
limited to tuned caoeadoat To the beat of dm 
authors* knowledge, the forced reopona’.!i of both 
tuned and mistimed cascades using multi*>degree 
freedom models has not been presented in the pub** 
llshcd literature. Thus, the problem considered In 
this paper is a logical extension to the present 
state of the iiterature on flutter and forced re- 
sponse of blades with mistiming# 

The mnUiematlcal model considered heroin is of 
the dlocreclRod, lumped parameter type, utilizing 
discrete nuisoos, mass moment of Inertia, and linear 
and rotational springs to represent the Indivlduol 
blades. The unsteady aemlynamie loads were calcu- 
lated by using liSiltehead*s(12) incompressible flow 
cascade theory. Thus, the model considered is 
simple enough to be used for extensive parametric 
studies* At the same time, it is adequate to repre* 
sent the basic dynamic characteristics of a mistimed 
cascade, to provide guidance In refining both the 
aerodynamic and structural models, ami to check the 
results obtained from finite element formulations, 
ouch as one presented in reference 13. Recently, 
the authors have extended the present work into the 
subsonic and supersonic flow regimes In refer- 
ence 14. 


n*,_jnAeor y 

Xn general, the components which comprise a 
bladed-dlsk system have complex geometries. The 
analysis of tills complex system, as stated earlier, 
Is further complicated by blade mistiming* To 
accomplish Che stated objectives of the paper, it 
Is necessary to develop a model vdilcli simplifies 
the analysis, yec maintains the basic dynamic char- 
acteristics. For this reason onlv two degrees of 
freedom (one bending ami one torsi ni) lor each 
blade are considered in this paper. However, the 
authors have plans to add additional blade degrees 
of freedom and disk flexibilities in addition to 
other refinements to the present model. The gen- 
eral motion of a mistimed cascade is assumed to be 
a combination of all possible mvUlons of the asso- 
ciated tuned '*ascade. It will, therefore, be 
instructive first to develop and understand the 
model of a tuned cascade. 


A . Tu nejJ Cascade J>lod e l 

The geometry of a tuned cascade model is shown 
in figure I, The disk is assumed to be rigid and 
the bladed assemblv is modeled as infinite two- 
dimensional cascade of airfoils in a uniform up- 
stream flow with a velocity V as illustrated in 
figure 1. The effects of wakes shed from upstream 
obstructions are included. The wakes considered 
are limited to sinusoidal distortions re pro sent ed by 
VO r t i c 1 C y p e r t u rba t i on , so that t he y a re convec ted 
downstream at the flow velocity V. The amplitude 
of the wakes is specified by the velocity which the 
wakes would induce at the position of the mid- 
chord point of the reference blade as indicated in 


figure 1# The motion of the alrfotlo In each mode 
of cho tuned cascadfi In nooumed to be simple 
harmonic with a constant pb^ioo angle between 
adjacent bindea. Aloe, this Intcrblade phase angle 
is reotiictod by bam *a (15) assumption to the H 
disc rate values gj. « 2nr/N where r » 0, I, 2 
N-1. Conoequently, there are N modeo for the 
cascade with each b^ade having the same amplitude, 
the motion of a tuned cascade In rth mode Involving 
bending and torsion coupling can be represented In 
the form of a traveling wave ns shovm in figure !♦ 
The motion of the sth blade when the cascade vi- 
brates in the rth mode is Indicated In figure 2. 

For a tuned oystem, the modes with different inter- 
blade phase angles are uncoupled and hence one can 
write 


l**ari 

Furthermore, it is adequate to analyze the motion 
of a single blade in each of the Interbladc phose 
angle modes separately* Hence, the number of de- 
greofl of freedom of a tuned cascade for the present 
case is reduced to two for each value of 


B* fiistimod^mK^ade Model 


In a randomly mlo tuned caotMJe, the blades arc 
not Identical and can have dlff event response ampli- 
tudes, In addition, the phase angle between mlja- 
Lv(\t blades can vary. Because of the spatial por- 
Idicity, the general motion of a blade in a mis tuned 
cascade can be expressed as a combination oi the 
motions in all poanlblo interblade phase angle modes 
ol the corresponding tuned cascade. Consequently, 
the motion of the sth blade can be written in the 
traveling wave fonn 


N-1 


ritai i.t V p.„\ 

kr “AU 

r^^O 


i (u)t+B,.o) 


( 2 ) 


The qumultlej haj.*s and a.„,*s were called as 
Che 'aerodynamic modes' in reference 4. For a cas- 
cade with N mintuned blades, equation (2) can be 
generalized as 


{XI 0^“*^ " [E](V) 0^“"^ (3) 
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It should be rera/itked that the motion of a 
blade In a caocado with and without mistuning can 
alDo be exproaoed i^ a otandlng wove form* Since 
the Gtability and tcoponoc are independent of the 
form uacd, cons' xderation of cither of theac forms 
is adequate to describe the motion. In moot of the 
published literature on flutter analysis of bladed- 
disk asoemblloo, the traveling wave form is prefer- 
red, This is in contrast to the conventional 
flutter onalyoto of a fixed wing aircraft in which 
the standing w.%ve form is generally used. The 
motion of the sth blade of a miatuned rotor in the 
standing wave form may be written os 



(5) 

where 

Oq r 0rs/r (6) 

It is of interest to note that the special coses of 
the form given by equation (5) are known as 'multi- 
blade coordinate transformations' in the literature 
dealing with the helicopter and prop-rotor aero- 
clast Iclty, This transformation has certain advan- 
tages in the analyses » particularly when the coupl- 
ing between the rotor and supporting structure are 
involved! and when the coefficients la the equa- 
tions of motion are periodic in time. Hence, the 
standing wave representation may be expected to 
hove slmildr advantages if the analysis includes 
bearing motion, whirling motion of shaft, stand 
motion, etc. An interesting discussion on both 
traveling and standing wave representations for 
tuned bladed disk assemblii s is presented in ref- 
erence 16. 


C. Structural Model 


The structural model of the sth blade of a 
raistuned cascade is illustrated in figure 2. Each 
airfoil is suspended by bending and torsional 
springs, and , respectively. The airfoil 


in aonumed to be rigid in tho eUordwine direatlortf 
and this motion in neglocted. The clastic coupling 
between bending and torolon dun to pi:tttwi 0 t» shrouds, 
and rotation of the rotor in modeled through the 
offset dlstanee (bx^^) between the center of gravity 
and clastic axis* The centrifugal ntlffealng ef- 
fects due to rotation arc included in the bonding 
and torsional spring constants. The elastic and 
dynamic properties of the blades arc represented by 
their respcGtivc values at the three-quarters sta- 
tion of blade span. This model may be viewed os a 
logical extension of the so-called 'typical section 
wing' used in fixed wing acroeloocicity. 


The unsteady aerodvaamic loads were calculated 
by using Whitehead's^ ^2) cascade theory in the In- 
compressible unsteady flow* This theory is an 
extension of two-dimensional unsteady airfoil 
thcc*’y of Theoderson to account for cascade offocts* 
The 0 -feet of airfoil thickness, camber, and t>* 4y 
state angle of attack are neglected. As meats .^ed 
earlier, the effects of wakes from a periodic 
obstruction upstream arc included in the form of a 
vort iclty perturbation# It should be noted that 
the direction of the velocity, as shown In figure 1, 
Induced by the wakes is opposite to tliat of refer- 
ence 12# In view of the basic objectives of this 
paper, it is felt tliat this incompressible theory 
is adequate. However, the comprcsnlbllity effects 
will bo included by using Smith' theory in the 
subsonic flow regime, and Adamezyk ond 
Goldstein’ theory In the supersonic flow 
regime# These results will be reported in 
refer ence 14, 


E, Equations of Motion 

A simple application of Lagrange's equation to 
the mathematical model of the sth blade in figure 2 
leads to the following coupled bonding- tors ion 
equations 


— 



IDg 

®“0 
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v.dt'^ 


(1 + 2i5hg)nigUihg 0 


a + 2lCag)l(,g«'„gj 



(7) 


Structural damping is added to the equations of 
motion by multiplying the uncoupled stiffness coef- 
ficients in the bending and torsion by (1 + 2itj^ ) 
and (1 + 2ii;Qj^), respectively# The critical damping 



ratiosi $1,^ and arc related to the logari- 
thmic decremontOj ond by the rclatioa 
d|, » 2n5jj and <? aerodynamic 

lorces due to motion^ ropreaented by the auper- 
ocripc Ml and due to cxcltotlon from Binuooldal 
wakes I represented by the superocrlpt w, aro ok- 
prosaed In terms of nondimenslonal coefficients os 
follows? 
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The coefficients (Cpq)y^i (Cp^)pi arc 

colculoted by the unsteady cascade airfoil theory 
of reference 12 for given values of k, b, s/c, C* 
ond a («2n - 1) » The quantity Wr is the ampli- 
tude of the velocity of the sinusoidal wake in the 
rth mode. Nondlmensionalizing equation (7), ex- 
tending the resultant equation to all the blades 
(s “ 0, 1, 2 ... N-1) * and using equation (3), the 
equations for all the blodes of a randomly mlatuned 
cascade can be simplified as 



Dpi - [n-OfY) “ 

where 


[PI ■= 

jjE]"^[D][E) + [E] 
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Ilia Solution 


The aeroelastic stability of the cascade is 
determined by eigenvalues, Y^o> of the matrix [P]* 
The relation between the frequency w and y le 


Iw/wq « i/ V y “ p* ± i^ (12) 

Flutter occurs when IF > 0. 


For the given values of the number of blades, 
and hence the allowable Interblade phase angles, 
the gap to chord ratio, the stagger angle, the 
elastic axis position, and the structural param- 
eters, the eigenvalues of the matrix [Pj are calcu- 
lated for a range of values of k. Denoting the 
values of k and v at which TJ » 0 as kp and 
■vp, respectively, the nondimensional flutter speed 
can be written os 


Vp/bu)Q “ >^p/kp (13) 

The aeroelastic response of the blades Induced 
by wakes is calculated from equation (10) and is 


(Y) ■=■ - [(P) - [I]Y]"^tE)”^tG)[E](AIl) (14) 


The amplitude of each blade is obtained by substi- 
tuting equation (1^) into equation (3). 


IV. Results and Discussion 


Computer Program and Verification 

A digltol computer program was written to cal- 
culate the flutter stability boundaries and the 
blade response of a randomly mis tuned rotor. In 
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thio prograra^ it in poooiblo to conoldct my typo 
of miatuning oweh m bjado«to«blailo vapiatlono ot 
tho unooupXctl bending and tornional froquenolcot 
damping ratiofif maoo ratloSf olaptic axta and 
confcot of gravity poaitionQ> and m on» Thio pro- 
grata in operational on the NASA bewia Kopoarch 
Center IdH 370/3033» tuned and mintuned uu*» 
coupled bonding and uncoupled toroion caoop, in 
addition to tho tuned coupled bending-toroion ease, 
can bo treated no opocial caoea of thin prograna. 
Thin program wao cheeked for tho following opocinl 
caaopt 

It The corrcctnoao of the aoiodynamie cooffi- 
cientn wap chocked by comparlpon of tho prooent 
rooulto to Whifiohoad’o publiohcd ropulto in rofor-* 
crtco 12 over a wide range of cascade paramotorp. 
Thio comparioon ohowp excellent agreement* 

2* To check tho corrcccncoo of ths uncoupled 
torsional eigcnvaluop calculated by tho present 
program, a I2**bladcd rotor doscribod by Whitehead 
in rofornco 3 was considered. Tho root locus of 
this rotor with mistuned blades was ptesontod by 
Whitehead In figure 18 of reference 3» A compari- 
oon of the results shows oxc client agreement except 
for the case of zero interblade plmse angle. This 
difference appears to be a plotting error in refer- 
ence 3 because the present rooults arc also in 
agreement with those of reference 7 for all Inter- 
blade phase angles. 

3. To chock tho correctness of the program in 
calculating tho coupled bending-torn ion flutter 
speed of a tuned rotor, the present results for a 
few puleoted cases were compared to the corresp nd- 
ing ones In reference 11. This comparison shows 
excellent agreement. 

4. Forced response results are very limited 
in the published literature. However, an upper- 
bound for uncoupled bending or torsional response 
was given in reference A. Tho reoultG of the 
present program are within that bound. 


B . Aoroolnstic Stability 

A compressor stage representative of a forward 
stage of an advanced axial flow compressor was 
chosen for conducting parametric studies. This 
stage, known as the NASA Test Rotor 12, is showi\ in 
figure 3. Tho required parameters of this stage 
were calculated from the data given in reference 20 
and arc listed in table 1. The blade bending to 
torsion frequency ratio for this rotor is 0.357, 
and the elastic axis and C.G. position arc at 50 
percent chord. As a result, the coupling between 
bending and torsion is very weak and tho flutter 
mode is dominated by torsional motion. Hence, the 
results for the predominantly bending modes for 
some cases will not be presented. However, to con- 
duct parametric studies the bending to torsion fre- 
quency ratio and eJastic axis position are varied. 
For these cases the results for prcdoroinently bend- 
ing modes ore included. In some selected cases 
only uncoupled torsional motion is considered. 

A comparioon of the system eigenvalues of both 
tuned and mistuned cascades is useful to understand 
the mis tuning effects. For example, figure 4 is 
such a plot for a special case in which only un- 
coupled torsional motion of the blades is consid- 
ered. The first type of mis tuning considered is 


the one iti wliich the odd and evm p'^mbered bladee 
have different torsional frequenciee> Thin la known 
an flltemate blade miefeuning. For oxamplsi In the 
caoe of one percent miatuning, the frequency ratio 
If 005 lor all the even blades and is 

0.995 for all tho odd bladeo# Tho roforonco fre- 
quency Id equal to the arltluaocie mean of the 
uncoupled torsional frequencies of all the bladoot 
Because of the symmetry of this type of mlstuning 
tho mode couples with the n) mode only* 
Figure 4 may be viewed as a root locus for tho 
tuned cascade with % as the parameter* When tho 
blades are mistuned, this doseription is not com- 
pletely appropriate because each mode contains all 
possible intorblado phase angles, Due to tho in- 
herent symmetry in alternate blade mlstuning, each 
mistuned mode contains only two interblade phase 
angle modes. However, one can view this plot as a 
root locus with a predominant interblade phase 
angle as a parameter. Several intercepting observa- 
uions follow from figure 4. Even me percent of 
mistiming significantly affected cho system aigon- 
values and stabilized an unstable tuned cascade. 

As the level of mlstuning is increased, the hori- 
zontal width of the root locus is decreased , Tnls 
amounts to eaylng that tho effettivo damping of 
some modes Is increased while that of others is 
decreased. This behavior will have an iJifluence 
on forced reoponoe which will b« »H«<*'^,osed later. 
Whoa the mlstuning level is 1.5 percent, the root 
locus is split into high and low frequency groups. 

In the high frequency group all the interblade 
phase angles between 0 and 90 and between 276.4 
ond 360 are predominant,* in the low frequency 
group those batwf»cn 96.4*^ ond 270^ are predominant. 
.As the level i ' mistuning is increased beyond 1.5 
percent, che frequency separation of these groups 
becomes larger end the area enclosed by each group 
decreoseo. 

."igure 5 shows the eigenvector corresponding 
to the leost stable point of the 1.5 percent mis- 
tuuing curve on figure 4. The predominant Intor- 
blade phase angle at this point is 51.42® (r « 8) 
and is coupled with 231.42® mode (r « 36). As can 
be seen, all the even blades have the same amplitude 
and all the odd blodes have the same amplitude, as 
expected from the symmetry of this kind of mis tun- 
ing. How war, the amplitude of the odd blades is 
68 percent of that of oven blades. 

Figure 6 Illustrates the variation of uncou- 
pled flutter speed with. the level of alternating 
blade mlstuning with and without damping. The 
value of tho damping ratio used is 0.2 percent and 
is the same for all the blades. The results indi- 
cate that tho mlstuning has a substantial effect on 
flutter speed. For the undamped case, the flutter 
speed increase? raonotonically with Increase in mis- 
tunlng level. However, when the mlstuning level is 
5 percent and m-jra, the additional benefit is 
modest. For tht damped case a similar variation in 
flutter speed is noticed, except when the mlstuning 
is between 0-1 percent. In this range, the damping 
effect is more pronounced. It should bo noted that 
these observations, particularly the effects of rois- 
tunlng, are in agreement with the qualitative con- 
clusions analytically reached in reference 4 and 
tho experimental results reported in reference 21, 

A second, mor<^ common type of mlstuning was 
analyzed. Blade torsional freruencies were ran- 
domly chosen from a normally distributed population 
with a mean of 1 and a standard deviation 


6 



of 0*005» fho resulting blade frequencieo are 
shovm f.> figure 7. h coaparison of both the tuned 
and Rlti^imod oigenvaluos ia ohovfn in liguro 8» Tho 
tuned oigenvolaeo are the name as those presented 
in rigure 4, As can ho soeni there la stabiliaing 
effect on t!m system, hut is not ^uite as ocrong as 
that produced by one percent alternating miotunliig 
shown in flaure A, The oiBonvactor for the least 
stable miotuned mode from figure B ia ahoi/n in 
figure 9* The eigenvector consiats prcdominanrly 
of 6A*a9^ Interhlade phase angle mode but also ban 
slgnlficont participation from the S7»86^^ 81*57®) 
19.29® and 77*14® modes. Alao, each blade has a 
different amplitude and the interblade phase angle 
varies considerably. It can be seen that this 
flutter mode can be viewed as a localised phenom- 
ena since only blades numbered 42 t* rough 51 have 
large relative amplitudes. This type of behavior 
lias boon observed in actual engine tests. 

the effect of claotio anis position on flutter 
speed is of interest, both uncoupled torslomil and 
coupled bendlng-torolon flutter analyses were 
formed with and without damping for three elastic 
axis positions with the center of gravity at mid- 
chord. Ao was noticed in teference 11, very weak 
irtotabllitleo appeared in some canefi which were 
eliminated by the addition of a omall amount of 
Dtructural damping* Consenuently, they are of 
little practical interest, and, hence, only the 
results with structural damping are presented in 
figure 10. Note, that for uncoupled torsional 
flutter the worst location of the elastic axis in 
at the midchord point (a « 0), Thin is In contrast 
to the present results without damping av'd to those 
In references 1 and ll in xdiich a severe drop in 
flutter speed is noticed for the 75 percent 
elastic axis position. For counled bending- tors ion 
flutter, the effect of elastic axis position de- 
pends or, (jjjig/u'o* if * 1. the best loca- 

tion (of the three locations investigated heroin) is 
the 75 poreont chord point; *f 

boot location io the ^5 percent chonl point. When 
the elastic axis in oft the midchord, the effect of 
bending-torsion coupling on flutter speed is slg- 
nlf leant. This observation in in agreement with 
that in references 5 and 11. These results suggest 
that the tailoring of the elastic axis position con 
be used ns a passive control to increase flutter 
speed ao is done in fixed wing aeroelastlcity. 

Figure 11 shown the effects of both alternat- 
ing blade mistuning and damping on coupled bending- 
torsion flutter speed. As can be seen, both mls- 
tunlug and damping have beneficial effects. How- 
ever, the level of benefit depends on 
should be notfui that the adverse effect of coupling 
between bending and torsion i., tuned cascade flut- 
ter speed when ^ reduced by the bene- 

ficial effect of small mistuning and/or damping. 

The effects of alternating blade mistuning, 
damping, and elastic axis position on the coupled 
bending torsion flutter speed arc illustrated in 
figure 12* The effects of mistuning and damping 
are similar to those discussed in figure 6j the 
effectfl of the elastic axis position arc similar to 
those dlDcunned In figure 10. For all practical 
purposes, the curves for the elastic axis at raid- 
chord are the same an those In figure f> because the 
ratio is small. 


The results presented thus far only show the 
effects of miatunlno on aoroelasti: ofcabiUty, 

These renal to suggest that the utilisation of mis« 
tuning and/or tailoring of the elastic axis position 
as passive controls to increano flutter speed are 
feasible* the next stop is to oxatjine the offeet of 
mistuning on forced reoponsc# 


In tho present formulatloni it is poasihlo to 
conoii’or on excitation function consisting of all 
harmonico of rotational speed of th^. rotor wTioh 
range up to r » N-l» In engine acroclastic term- 
inology, the harmonic number r to known as the 
‘englno order* of ,hc excitation. The coefficients 
%ir ^w'.»r enuatlono 8(a) and (b) represent 
the forcing funetienc in the bending and torsion 
equations, respectively* To understand the nature 
oi the resoonoe, excitation in only one Immonie at 
a time will be conoidered* Thin results in no loss 
of genorality obtee the principle of nuporposition 
holds. If the r » II harmonic is considered, then 
the column nuitrices CAD^I, CAD|h {ADn«il are 
gero except I ADgl In equation til) . This corres- 
ponds to the ease in which there are R oyrametri*’ 
rally spaced obstructions located upstream from the 
blades and the circumferential wako distribution io 
perfectly sinusoidal. For practical applications, 
tho forcing frequency io thus equal t'j H times 
the rotational speed, 

The aeroeiastic responsa results presented 
herein are or two values of U, 11 and 39, at a 
fixed reduced frequency almsen such that the cascade 
is aeroelancicaXly stable in all modes. These 
values for H were packed because the oerodynamia 
damping of tlie tuned nyotem in tho r « 11 mode is 
relatively low wticreao that in the r « 39 mode is 
relatively high. The forcing frequency range in- 
vestigated is limited to a small range around tive 
uncoupled torsional frequency, 

If the blades are. tuned, the response will be 
entirely in the t ^ l\ mode, and all the blades 
have equal umplltude. The amplitude of response of 

J 

any blade in 1 which io a function of w/oso* 

let the torsional amplitude of resonance of each 
blade of the ;uned rotor be If the blades 

are now rarJomly mis tuned, there will be a response 
in all t\*i modes (enumerated by r) and the ampli- 
tude of response of the sth blade in 



Figures 13(a) and (b) with R « 11, and figures 
13(c) and (d) with R «» 39 show tho variation of 
®s/%,id ^®*^h the tuned and one percent alter- 

nating blade mis tuned cascades. Note that figutes 
13(b) and (d) are a repetition of figures 13(a) ond 
(c) , respectively, with 0.2 percent structural 
damping. Note that the torsional amplitude 
of the tuned eaoeade depends on the level of damping 
and R. The bonding ampUtudoo ai\^ not ohuwn because 
they are very small In the range of the excitation 
frcQuency shown herein. For the alternating blade 
mistuning, only two r modes ate coupled. Tho 
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flsplifcutjiQ behavior is siaiiar to fcJiai; Qho*»m In 
uro ^ in ohich t\u acpXitu^ea of the oM amJ oven 
biadea are different* In all these cases» the 
oin^la resonance peak of the timed cascade is re» 
placed W t*Hn reoonanee peaks for the alternatlnB 
mistuned cascade* It is mm tJiat the effect of 
nlotuning on forced response depends on the engine** 
order of the lorcioB function* For eHajuploi the 
lalotoning lias a beneficial effect (fiB* l3Cs)) on 
torsional response for II » ll, but has an adveroo 
effect (fin* l'J(c)) for R » 39. this lo in con- 
trast to the coimTrOn belief that the mintunlng 
always has an adverse effect on forced response* 
Thus, this result provides an added incentive for 
purauing the use of mistiming m a passive control* 
The maHimum docr^ase in amplitude with mistiming 
for R ^ 11 is appi’OKlmatoly 89 percent (fig* 

13(a)) without damping and is approKimately 2S per- 
cent (fig* 11(b)) with damplngf The maximum in- 
crease in amplitude v;lth mistuning for R * 39 ie 
approKinmtely IlO percent without damping and is 
approKimately iO percent with damping* As expected, 
the nr ’vetural damping has a significant effect on 
forced response. Although not shown, the decreaoe 
in tuned issonance response is 01 percent for 
R « ll and 37 percent for R « 39* 

The randomly mistuned cascade described in fig- 
urn 7 was analy;5ed for forced response with k« 1.2. 
The rcnulto are ohoim in figures l^4(o) for R U 
and 14(b) for U 39. As in figure 9, each blade 
has a different amplleudo and the Interbladc phase 
angle varies. It is seen that the single resonance 
peak for the tuned case is clumged inf* . *ultlpie 
peaks, and different blades peak at d"l' /rent forc- 
ing fretiuoncloo* As in the alternating mlo tuning 
case, this type of miotuning has a beneficial effect 
on the R » 11 excitation. For the R 39 exci- 
tation, the sharp resonance peak of the tuned oyo- 
tem Is eliminated. However, for moot of the values 
of the forcing frequency, the mistuned response is 
higher than the tuned responoe. 


V. Conclusions 

This investigation was conducted in on attempt 
CO improve the basic understanding of the effects 
of mis tuning on aeroclaotlc stability and rcoponGo 
and rhen to explore the feanlbility of using mla- 
tuning as a passive control to increase flutter 
speed mid to minimize response. The following con- 
clusions are reached on the basis of the limited 
results obtained by using incomproooible unsteady 
cascade aerodynamic thooty for two typos of mis- 
tunings 

1. In general, the mistunlng has a beneficial 
effect on the coupled bending-torsion flutter speed. 
The flutter speed increases monotonically with an 
increase in alternating blade mis tuning level. How- 
ever, when the mlstuning level is above about 5 
percent and more, the additional benefit ia modest. 

2. The inherent random mistiming which exists 
in real fan, compressor, and turbine stages has a 
olgnif icontly beneficial effect on flutter. This 
observation Is qualitatively in agreement with the 
experimental results published in the literature. 

3. As expected, the effOLt of structural damp- 
ing on flutter speed is stabilizing. However, In 
the presence of mlrtunlng this effect is not as 
significant as in :he tuned case. 


4» The uao of uncoupled toraional flutJtcr 
snalyoio to deduce the effect of olantic axis posi- 
tion was found to bo unrollablo feeemmo the 
ing between bending and torolon, structural damp- 
ing, and mlstuning can change the rosulto 
oignlflcani^iy* 

St Mlsfuning have either a bmeficial or 
an adverse es fvet on forced reoponno, depending on 
the engine ordu% of excitation* 

6* ’“istunlng introduced multiple rusonant 
peaks for a given euoino order excitation. 
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Table 1 Faramctcro of NABA Toot Kotor 12 


N 

55 

fi/c 

0.514 

>'o 

258.5 

a 

0 (varied in oome caoeo) 

^to 

0 (varied In some caoca) 


J 0 5774 (.; ° 0) 


1 0.7638 (0 « -0.5 nnd 0.5) 

5 

J4.4° 

^jg/wajj (tuned) 

0.357 (varied in oome caBoa) 
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Figure 1, - Geometry of a tuned cascade in rtii mode (note that the variable 
s represents the gap between chords as well as the blade number Index). 




Figure 3. • NASA test Rotor 12. 



Figure 4. - Effect of alternating blade mistuning on 
eigenvalues (torsion onlyl; a • 0, k • 0,642, 
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Figure 5. - Effect of alternating blade mistuning on the eigenvector of 
the least stable eigenvalue from Figure 4\vith 1 % mistuning. 



Figure 6, - Variation of uncoupled torsional flutter speed 
with aiternating blade mistuning both with and without 
damping; a = 0. 




BLADE NUMBER, s 

Figure 7. - Individual blade torsional frequencies for 
the randomly mistuned case. 



REAL PART OF EIGENVALUE, p 

Figure 8. - Effect of random mistuning on eigenvalues 
(torsion only); a - 0, k = 0. 642; - 0. 










TORSIONAL AMPLITUDE RATIOS. Oj/a^ jjj 



(b)R = ll. 0.002. 

Figure 13. - Effect of blade alternating mistuning on 
coupied bending torsional response-, a • 0, 

W|^ /Wq° 0. 357, k - 1.2, c.g, atmidchord. 
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